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Let H(t) = (1 — t/T)Ho + (t/T)H\ , t £ [0,T], be the Hamiltonian governing an adiabatic 
quantum algorithm, where Hq is diagonal in the Hadamard basis and Hi is diagonal in 
the computational basis. We prove that Ho and Hi must each have at least two large 
' mutually-orthogonal eigenspaces if the algorithm's running time is to be subexponential 

S-H ' in the number of qubits. We also reproduce the optimality proof of Farhi and Gutmann's 

search algorithm in the context of this adiabatic scheme; because we only consider initial 
Hamiltonians that are diagonal in the Hadamard basis, our result is slightly stronger 
than the original. 

Keywords: adiabatic; quantum; bounds. 



1. Introduction 

Farhi et ali^ proposed a continuous-time quantum algorithm for solving NP- 
complete combinatorial search problems, invoking the adiabatic approximation!^ It 
has been shown that Farhi et al.'s original proposal for the 3SAT problem has expo- 
nential worst-case complexity!^ However, the average-case perf ormance of (possible 
modifications of) the original proposal is still an open question N^ I ^ A more gen- 
eral notion of adiabatic quantum computation (n ot consid ered here) is universal.^ 
This worljfl supplements the similar recent worlJ2liniIl] on analytic lower bounds 
on the runtime of certain classes of simple adiabatic quantum algorithms. Our as- 

a previously appearing in Ref. 1161 

1 



I 



I 



June 2, 2008 15:4 WSPC/INSTRUCTION FILE IoannouMosca'IJQI 



2 LAWRENCE M. IOANNOU and MICHELE MOSCA 

sumptions are weaker: for our general bound, we do not assume that the initial 
Hamiltonian is a one-dimensional projector. Thus, the dimension of the second- 
largest eigenspace of the Hamiltonians emerges as a factor in the algorithms' com- 
plexity. 

Let the function / : {0, 1}" — > {0, 1} such that 

10 it z = w 

encode a general search problem (with unique solution w) for which we seek a 
quantum algorithm. The uniform-amplitude superposition state, 

W = -m E l*>> ( 2 ) 

VJV ze{o,i}™ 

will be the starting state (at t — 0) for all the algorithms we consider. Define the 
Walsh- Hadamard transform (or just Hadamard transform) 

W = ^m E E (- i r v \x)(y\- (3) 

^ x£{o,i}« j/e{o,i}» 

Let 

|5)=W|«), for all z G {0, 1}". 

We define the Hadamard basis to be {|z) : z £ {0, 1}"}. 

Recall from Ref. HU and Ref. [13] that the continuous-time version of the general 
search problem, defined in Ref. 14, may be solved adiabatically with initial and 
final Hamiltonians 

H u = ^2 h scllIch (z)\z) (z\ (4) 

2G{0,1}" 

and 

H w = f{z)\z){z\, (5) 

zG{0,l}" 

respectively, where h scal - ch (z) £ {0, 1} and h scarch (z) = & z = 0". 

If / is the truth function corresponding to a Boolean formula (also denoted 
/) which is an instance of the 3SAT problenff^, then we can define the function 
v : {0, 1}« — ► {0, 1, . . .}, 

v(z) = the number of clauses of / violated by assignment z. (6) 

Recall that Farhi et al.'s original adiabatic algorithm for 3SAT uses initial and final 
Hamiltonians 

ze{o,i} n 
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and 

H F = «(*)I*X*I> (8) 

ze{o,i}™ 

respectively, where h(z) = J2i- Z =i ^* an< ^ ^» ^ s * ne num ber of clauses of / that 
contain variable 

2. Results 

It is convenient to cast Farhi et a/.'s adiabatic 3SAT algorithm in a more general 
notation. Let N — 2™. Let p(n) and q(n) be two nonnegative integer functions of 

n. Let V = {ij-}j=o,i,..., p (n) and 2 = {<2fc}fc=o,i,. ..,?(") be two P artitions of {0, i}™, 
that is, the following two unions are disjoint: 

q(n) p(n) 

|J Q fc = {0,l}»= (JP,, (9) 

fe=0 j=0 

Define the projectors 

g fc ^^|z)(z| and £,-= £|*><z| (10) 

onto the eigenspaces of the two Hamiltonians H and Hi , 

<?(«) P{n) 

H Q = Y,F k Qk and H^^EjPj, (11) 

where the two sequences of energy eigenvalues (-Ej)j=o,i,..., P (n) and (-Ffc)fc=o,i,...,g(ra) 
are strictly increasing sequences of real numbers. Assume further that F — = Eq 
and that the starting state \u) is a ground state of H , i.e. 0" € Qo- We will also 
write 

H = F ( z )\ z )( z \ and Ht = ^ E ( z )\ z )( z \ ( 12 ) 

ze{o,i}™ zG{0,l}" 

for eigenvalue functions E(z) and F(^) which are defined by the partitions V and 
Q and the sequences (Ej) and (-Ffc). 

Let s(f) be a smooth, nondecreasing, real function of time t such that 

s:[0,r]— »[0,1], (13) 

with s(0) = and s(T) = 1, for some final time T. It is clear that 

H gen (t)= (1- a(t))H + a(t)H u 0<t<T, (14) 

is a Hamiltonian having suitable form for adiabatic quantum algorithms solving 
problems for which Pq contains the solutions. 
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2.1. A general bound 

Let |'0 g en(i)) be the state evolving under H gen (t), that is, we define the initial value 
problem 

J^genW) = -lffgen(*)|Vgen(i)>, <t <T (15) 
IV'gen(O)) - \U). (16) 

Let K be the index of the largest set Qk- 

|Q«-|>|Qfc| for aUfc = 0,l,...,g(n). (17) 

Note 

/ \ 

Hgen(t) = (1 - «(*)) Uw/-5Z( i ^- J? *)Qfc +«(*)S^)kX«l- (18) 

\ fe=0 / 2 

Noting that 

m = w\ Z ) = -Lj2Y,(- 1 ) xmv \ x )(y\\ z ) = -t^Ew*^ ( 19 ) 

we have, dropping the time-dependence notation, 

#gc„ = (i-a) (f k i - ^ E E ( eW - E (-i)" (xew) ) k>(yi 

+ aY,E(z)\z)(z\. 



Suppose |V'gcn(t)) is expressed in the computational basis as IV'gen(i)) = 
^^J^ 1 7 z (t)|z). Clearly, the goal of the unitary evolution is to make 
T, w 'eP IWIV'genCT)}! 2 = T, w 'ePo l7^'( T )| 2 of constant (or just l/poly(n)) order 
so that we have a good probability of discovering one of the solutions w' by per- 
forming a measurement in the computational basis on the final state \ip gen (T)}. The 
Schrodinger equation implies 



J^gcn) = ~*E (J2 ( X \ H Scn\yhy) \x) . 

x \ y / 
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For each solution w' G Pq, 



9 I |2 



2Re(( W '| J^gcn) • Wgcn}*) 



2^{«'| J ff gen |y}Im ( 7l0 , 7 ;) 



1 9(n) / \ 

-2(! - ^ E E( F * - ^) E c- 1 )"^ ey) 

y ^w' k=o \zeQk J 



Im (7^7* 



q(n) 



E (-ir (w ' ffia) ] im ( 7 »o 



= -2(i--)^EE(^- F *) 

where the last two lines hold because Im^Tu/T^/) = 0. If the strings z G Qk were 
somehow "random enough", we would expect the sums (j2 ze q k (— l) z *( w ® yS >^j to 
be close to 0. This would suggest that the derivative is small for an average problem 
instance. Without going into statistics about typical problem instances, we can only 
bound the size of these sums by \Qk\ - Using this bound for every w' G Pq and using 
Im( 7l u/7*) < y/N (which follows from applying the Cauchy-Schwarz inequality 
to the absolute value of the dot-product of the vector of all l's and the vector of 
Im(7^'7*) values) we get 



y - 

^ dt 

w'eP 



In 



< 2 



\P 



■ maxffc • dim (span{|z) : z G 
TV ^ 



IK 



(20) 



Integrating from t = to t = T as before, in order to lower-bound by a constant c 
the total probability of finding a solution, we require 



T > 



N-\P \/VN 



(21) 



2|P | ■ maxfe F k ■ dim (span{|z) : z G Qk}^) ' 

Assume that the number of solutions |Pq| is small. Thus if the dimension of the 
subspace orthogonal to the largest eigenspace of Hq is small, then the algorithm 
requires exponential time T. 



2.2. Searching with a more complex initial Hamiltonian 

We can use a similar argument to prove the analogous result for H\. In this section, 
we illustrate the main ingredients of this similar argument by giving it in the simpler 
case where the final Hamiltonian is a one-dimensional projector. Thus, as explained 
below, we get an interesting variant of Farhi and Gutmann's optimality result in 
Ref. HU We note that an argument similar to that in Ref. OH based on the time- 
reversibility of quantum mechanics, can be used to prove the following result; that 
is, the following result is implied by the "time-reversal" of the result in the previous 
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section (or the results in Ref. [9] or Ref. [TOl both of which consider one-dimensional- 
projector initial Hamiltonians). Our short proof below need not appeal to time- 
reversibility. 

Consider the following problem: "What is the minimum time T needed such 
that the Hamiltonian 

H B (t) = H gCD (t) | p(n):=li | Po |. =1 (22) 
= {l-s{t))^2F(z)\z)(z\ + s(t)E 1 (I-\w)(w\), 0<t<T (23) 

z 

(adiabatically) evolves the start state |?/>(0)) = \u) to a final state \ip (T)) that is 
close to |u>)?" Noting that 



x y 



we have 



\ x y 



(24) 



(25) 



Letting \ip(t)) = ^2 z /3g(t)\z), the Schrodinger equation implies that 

j^) = -iH B \^) = -xEE^i^iy)^!^)- 



x y 



Instead of bounding the rate at which probability amplitude flows into the solution 
state \w), we bound the rate at which it flows out of the start state \u). Noting that 
(u\ip(t)) = Po{t), we can derive that 



d_ 
dt 



2sE x 
N 



j2(-ir- v MPyPo) 



Using the Cauchy-Schwarz inequality we get the bound 

2Ei 



Noting that /3q(0) = 1 and using 



< 



N 



|/3o(T)| 2 -|/3 (0)| 5 



d\MW 
dt 



< t < T. 



dt > 



r 


d\Po(t)\ 2 


10 


dt 



dt 



gives 



IA)(T)| 2 >1- 



2E X T 



N 



The following inequality holds ff^ 

\(w\iP B (T))f + \(u\i; B (T))f<l + \(u\w)\. 



(26) 



(27) 



(28) 



(29) 
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This implies 

l( * (T)>|2 ^w + ^ (30) 

which implies that the algorithm requires exponentially large time T in order to 
compute w with sufficiently high probability. 

The result effectively bounds the power of any initial Hamiltonian that is di- 
agonal in the Hadamard basis, in the presence of a final Hamiltonian that merely 
encodes the binary function f(z) defining the generalized search problem (with 
unique solution w). We interpret this result in the context of Farhi and Gutmann's 
continuous-time search algorithm!^ Recall that their optimality result is that no 
Hamiltonian of the form Hp (t) + H w can solve the continuous-time search problem 
faster than H u + H w if it must work faster on most problem instances w £ {0, 1}™. 
It follows from our result that no Hjj(t) that is restricted to being diagonal in the 
Hadamard basis can outperform H u on even one problem instance. 

Conclusion 

Combining the techniques of Sections 12.11 and 12.21 we can show that, in our gener- 
alized adiabatic scheme, both Hq and Hi must have at least two large mutually- 
orthogonal eigenspaces if T is to be subexponential in n. 
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